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Abstract 

We systematically compute and discuss meson and muon polarized radiative 
decays. Doubly differential distributions in terms of momenta and helicities of the 
final lepton and photon are explicitly computed. The undergoing dynamics giving 
rise to lepton and photon polarizations is examined and analyzed in the soft and hard 
region of momenta. The particular configurations made by right-handed leptons 
with accompanying photons are investigated and interpreted as a manifestation of 
the axial anomaly. The photon polarization asymmetry is evaluated. Finiteness of 
polarized amplitudes against infrared and collinear singularities is shown to take 
place with mechanisms distinguishing between right handed and left handed final 
leptons. We propose a possible test using photon polarization to clarify a recently 
observed discrepancy in radiative meson decays. 



1 Introduction 



Radiative decays of light mesons and leptons have been widely studied both experimen- 
tally and theoretically. They represent an excellent source of information on the experi- 
mental side as well as a benchmark for theoretical speculations. Extensive comparisons 
have been carried on in the past between experiments and theoretical predictions for me- 
son radiative decays (see for example Ref.[l]). A while ago, radiative polarized leptonic 
decays of mesons [2,3] and muons [4,5] have also been considered. Recently special at- 
tention has been given to the role played by the final lepton mass mi in the threshold 
region of the decay and to the mi ^ limit concerning the helicity amplitudes for mesons 
[3] and leptons Refs.[4, 5]. The 0{a) radiative corrections generate an hehcity flip of the 
final lepton even in the zero mass limit [6] provided the lepton mass is kept from the 
beginning into account. Following the interpretation due to Dolgov and Zakharov [7] of 
the axial anomaly the final states with opposite helicity can be interpreted [3,8,9], as a 
manifestation of the axial anomaly giving rise to a peculiar mass-singularity cancellation 
for the right-handed polarized final lepton amplitudes. 

We consider in this work the case of polarized radiative decays of the pion and kaon 
meson and of the muon more extensively by taking into account polarizations of final 
lepton and photon degrees of freedom. Contrary to the previous case [2] , in meson decays 
we consider the polarization states of both lepton and photon final states. 

This approach, containing a complete description of the final momenta and helicities, 
may give further and more detailed information on the final state with respect to the 
inclusively polarized and unpolarized cases. Furthermore, the agreement with the more 
inclusive results previously obtained in the literature can be easily recovered by summing 
over the emitted final states polarizations. It is worth noticing that this approach allows 
to describe more closely the interplay between several peculiar features of the dynamics 
involved. As, for instance, to pinpoint the role played by angular momentum conservation 
and its connection with hard and soft photon momenta, and to consider the role played by 
the parity conservation in weak decays. All these aspects related to angular momentum 
dynamics may be effectively described in terms of the photon polarization asymmetry. 

Here we emphasize that the knowledge of the helicity amplitudes of the final leptons 
and photons, in addition to an explicit test of the angular momentum conservation, shows 
the relative rates of the partial helicity amplitudes. Indeed, in the total rate, difi^erent 
helicity amplitudes, depending on the range of momenta, enter with varying weights. 
Therefore, this behavior gives the opportunity to isolate peculiar polarized configurations 
in order to maximize or minimize them according to favorable intervals of momenta. As 
far as phenomenological applications are concerned, this may be, as will be discussed 
later, an effective way to compare theory and experiment on a new basis. The case of 
the photon polarization asymmetry, proposed in this work, allows, in this respect, a new 
approach to inspect interaction dynamics via a finite and universal quantity which is also 
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directly associated to parity violation. Moreover, the photon polarization asymmetry is 
very sensitive, in radiative meson decays, to the hadronic structure, allowing for a more 
precise determination of the electromagnetic form factors with respect to the one obtained 
so far. 

Some of the results achieved in this paper can be shortly listed: we explicitly calculate 
amplitudes and final distributions in terms of lepton and photon momenta at fixed final 
lepton and photon helicities. Double differential expressions in terms of lepton and photon 
momenta are also provided together with the partial helicity amplitudes for the meson 
and muon cases respectively. Moreover, we analyze how the cancellation pattern of mass 
singularities works on polarized processes. In the inclusive quantities this is a sensible test 
of the consistency of the results. Once inclusive distributions are obtained by integrating 
over final momenta, we observe the cancellation of all mass singularities both infrared 
and coUinear. In particular, a peculiar pattern of mass singularity cancellation is shown 
to take place, which differs for the left-handed helicity final lepton states with respect to 
the right-handed ones. The same behavior can be observed for the meson as well as for 
the muon case. 

Finally, we discuss a possible interpretation in terms of tensorial coupling of the results 
obtained recently at the PIBETA experiment [10] for the radiative pion decay in electron 
channel. It is argued that polarized radiative processes may constitute a sensible test to 
resolve the controversial issue of tensorial couphngs in radiative pion decay, allowing also 
for a sensitive test in the corresponding kaon decays as well. 

The paper is organized as follows: In Section 2 we consider the case of the meson 
polarized radiative decay. We discuss the contributing amplitudes and the underlying 
theoretical tools. We also define the gauge invariant set of matrix elements together with 
the definition of the Lorentz invariant quantities. Allowed and forbidden helicities con- 
figurations are here analyzed as well. In Section 3 the polarized radiative muon decay 
is discussed. In Section 4 we define distributions of branching ratios in the photon and 
electron energies and the photon polarization asymmetry. Numerical results for the dis- 
tributions of branching ratios and polarization asymmetries are provided in subsections 
4.1, 4.2, and 4.3 for the cases of pion, kaon, and muon decays respectively. Results for the 
polarized electron energy spectra are presented in Section 5. In Section 6 the dependence 
of the photon polarization asymmetry, induced by tensorial couplings, is discussed in the 
radiative pion and kaon decays. The peculiar pattern of mass singularities cancellation 
in polarized radiative decays is described in Section 7, while conclusions are presented 
in Section 8. General results and the corresponding formulae for the polarized radiative 
meson and muon decays are collected in Appendix A and B respectively. 
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Figure 1: Feynman diagrams for {tt^,K^) vil'^^ decay, wliere I = e, iJ, 

2 The polarized radiative meson decay 

We start this section with the calculation of the polarized amplitude for the process 

M^ip)^uiip,)l+ipi,\i)l{k,X,), (1) 

where M+ = 7r+ (K^) and / = e (/i) stand for pion (kaon) and electron (muon) respec- 
tively, with z/;=e,At ^hc Corresponding neutrinos. The four momenta p,Pu,Pi correspond 
to meson M, neutrino, and charged lepton, while A;,A^ indicate the charged lepton and 
photon helicity, respectively. The neutrino is assumed massless and therefore is a pure 
left-handed polarized state. The Feynman diagrams at tree-level for this process are 
shown in Fig.l, where the green bubble just indicates the Fermi interaction. The first 
two diagrams Figs.la-b correspond to the so-called inner bremsstrahlung (IB) diagrams, 
where the photon is emitted from external lines and the meson behaves as a point-like 
scalar particle. The third diagram Fig.lc is the so-called structure-dependent (SD) dia- 
gram, where the photon is emitted from an intermediate hadronic state and the matrix 
element will depend on the vectorial (V) and axial (A) meson form factors. The total 
amplitude for this process can be split in two gauge invariant contributions 

M(^i.^7) = M^^'^^^ + Mg^'^^^ , (2) 

where Mjb and Msd correspond to the IB and SD part of the amplitude. The IB ampli- 
tude is given by [11, 12] 



(3) 



where Jf, = 7°/ca, e^(fc, A^) stands for the photon polarization vector of momentum k and 
helicity A^, while u{p^) and v{pi, A;) are the bispinors of final neutrino and charged lepton 
respectively. Explanations of other symbols appearing above are in order. The is the 
Fermi constant, mi is the charged lepton mass, /m is the meson decay constant, where 
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/tt ~ 131 McV and fx — 161 McV, and Vug is the Cabibbo-Kobayashi-Maskawa matrix 
element corresponding to u ^ q = d and u ^ q = s quark transitions for pion and kaon 
decays respectively. 

The SD part of the amphtude contains vectorial (V) and axial {A) form factors, that 
clearly depend on the kind of initial meson, but not on the lepton final states. Indeed, 
they are connected to the matrix elements of the electromagnetic hadron current and 
the axial and vectorial weak currents A'^ and respectively, as 

{Y,AY^{p,k) ^ J dW^'^(0| TV,^{x){V{0),A{0)r \M^{p)) • (4) 
Using Lorentz covariance and electromagnetic gauge invariance, it follows that: 

rriM 

A.^fe .) ^ . .) A - - /„ + '-^^^) . (5) 

where r]'^" = diag(l, — 1, — 1, — 1) is the Minkowski metric, and e^^"'^ is the totally anti- 
symmetric tensor ^. Finally, the SD part of the amplitude is given by [11, 12] 

V2 ^ [ ruM V [p-k)j ruM J 

X [uiPuhi. (1 - 75) v{pi, Xi)] , (6) 

where ttim stands for the meson mass, while V and A are the meson vectorial and axial 
form factors respectively. Notice that both the terms Mjb and Msd are separately gauge 
invariant, as can be easily checked by making the substitution e*(A;, A-y) — > e*(A;, X^) + 
in Eqs.(3) and (6). 

Now we provide the corresponding expressions for the polarized amphtude in the 
center of mass (cm.) frame of the fermion pair (neutrino and charged lepton), namely 
Pi + Pu = 0. We choose a frame where the 3-momenta of neutrino and photon have the 
following components in polar coordinates 

— * 

Pt, — El, (sin ^ cos </?, sin ^ sin (/?, cos ^) , Pi — —Pu, k — Ej {0,0,1), (7) 

where E^ and E^ are the neutrino and photon energies respectively and 6, ip are the usual 
polar angles. For the photon polarization vectors we choose helicity eigenstates {e{k. A)), 
which in this frame are given by 

e^(A;,A^) = -^(0,l,a^,0) (8) 



^In our notation, the ^^^^"-0 is defined as e°^^^ = 1 and €0123 = — 1, when generic four- vectors are 
= {vq,v) and v^^ = {po, -v). 



4 



whose helicity eigenvalues correspond to A^, = —1 left-handed (L) and A^, = 1 right-handed 
(R) circular polarizations. Photon polarization vectors satisfy the transversality condition 
k^ei^{k,X^) = 0. Regarding the polarization vectors of fermions, it is convenient to use 
the solution of the Dirac equation for the particle (u) and antiparticle (v) bispinors in the 
momentum space [13]. In the standard basis ^ we have: 

where the 2-component spinors u!x{n) (with helicity A = ±1) are the eigenstates of the 
helicity operator {a ■ n) u!x{n) — Xu!x{n), and are the Pauli matrices. Here, n = p/\p\, 
where p is the 3-momentum and E — ^J\'p\'^ + m? is the corresponding energy. If p = 
|p| (sin ^cos </?, sin ^sin (/?, cos^), then in polar coordinates, uj\{n) can be expressed as 

cos ~ \ ( 6 sin — \ 

i£ • ' '^-i(^) = iv ' ] ■ (10) 

2 sm I y y 2 cos | y 

At this point it is convenient to introduce the following Lorentz invariant quantities 

_2p-k 2p-pi 2pi-k 

x= — —, y= — —, z = — —=y-l + x-ri (11) 
nij^ rrij^ rrij^ 

where in the meson rest frame, x and y are just proportional to the photon and charged 
lepton energies respectively and r; = mf/m\,j. Finally, after a straightforward algebra, 
the IB and SD contributions to the polarized amplitude in the fermion pair cm. frame 
are given by 

M^y-'^ = eGFmifMVu, - {Sx„-i (5x„-i + E^ R+ sin 9 
+ 5x„+i 5x„-i E^R- (1 - cos^)} e^^-<^ 

eGFfnljVuq ^^^"^^ Sx^,±i xj =f 5a, -i i?_sin6' 
± 5x„+iR+{cose±l)} e'^'''^, (12) 
where the structure dependent part is given by M^^^''^''^ = M^'^^'+''^ -|- M^'^^'^''^ and the 



TV/r(A;,A7) 



symbol R± = y Ey yy Ei + ^/Fi ± y Ei — y/nj , with Ei = Ei/niM and Ei is the energy 
of the final charged lepton. In this frame, the energies normalized to the meson mass are 
given by 

- X - 1 — X — ri - 1 — X + ri 

TP 77' '77' ' 



2^1^' ^ 2v/r^' ' 2v^ 



X 



cos 9 = (x-2)(l-x + n)+2a(l-x) 

x{\ — ri — x) 



^In the standard basis representation, 70 = Diag(l,— 1), and 7 = ( n I' 75 = ^ ^ 
where 1 = Diag(l, 1) and a are as usual the PauU matrices. 



-CT /' V 1 
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and 




= Vl-ri-x, R_^Jri^ . (14) 



Notice that, as expected from general arguments, the azimuthal angle (f factorizes in the 
overall phase of the amphtude. At this point it is important to stress that the SD terms 
in the amplitude, proportional toV + A and V — A, correspond to pure right-handed and 
left-handed photon polarizations respectively, while the IB one is a mixture of both. In 
particular, the terms proportional to pure left-handed photon polarizations in the Mib, 
come only from the tcnsorial structure in Eq.(3), namely from terms proportional to 
[uu -l- 75) vi], while scalar contributions of type [mj, (1 -|- 75) vi] do not select any 

specific photon polarization. We will return on this point in the following when anomalous 
tensorial coupling in radiative pion, and kaon decays will be discussed. 

By using Eqs.(12) and (13), it is now straightforward to evaluate the square modulus 
of the amplitude. Below we will provide its expression summed over the charged lepton 
polarizations, as a function of the photon helicities. After integrating over the phase 
space, we obtain for the photon polarized decay rate F^t^, the following result: 

f;;^=^ElM--f^.-(^.A). (15) 

Here itlm stands for the generic meson mass mM=n,K, and A = z/x. The Dalitz plot 
densities p^^ {x, A) for the polarized decay are Lorentz invariant functions, and are given 

by 

p^-'\x, A) = A, fi,(x, A) + Asn liV- Affj,{x, A) + A,^,{V - A)f}^,,{x, A) (16) 



p(+^)(x. A) = A, ff,{x, A) + Asn + A)'ff,{x, A) + A,^^{V + A)ff,,{x, A) , (17) 
where 

A - 9r,(l^\ A A -9r,I^A 

\mMj ruM 

Asn = ^G|mi,|K.r- (18) 

In the following, for later convenience, we will introduce the labels R and L corresponding 
to photon helicities A^ = 1 and A-y = —1 respectively. The functions f}"g^{x, A), fsh^{x, A), 
and /j^t(^) ^) given by 

/.'.(^,A) = l--^(i + r,(x-l)-^(l + x-n)) 
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ff,{x,X) = x'X((i-x){xX + ri)-ri) 

fj,{x,X) = x'{l-X){{x-l){ri + x{X-l))+ri) 

/L^^(x,A) = l^(x2^(i-x)(xA + n)-n) . (19) 

The function fj'gix, A) + f^^ix, A) coincides with the corresponding IB function fisix, A) 
for the unpolarizcd case provided in [11, 12, 14], as well as fsoix, A) = fs^^^x, A) + /j^(a;, A) 
and fiNrix, A) = f^j^^{x, A) + f^j.{x, A). More general results for the complete polarized 
radiative decay rate, including also the charged lepton helicity in the pion rest frame, are 
provided in appendix A. 

In order to obtain the differential branching ratio (BR) it is convenient to rewrite the 
term Asd in Eq.(18) as 

^.. = ^-7^(9^)'ro (20) 

47rri{l-ri) \ Jm J 

where Tq — r{M Ivi) is Born contribution to the total width of non radiative decay 
M — > Zf;, in particular 

ro(M+(p) ^ + /+) = (1 - r,f . (21) 

Then 

d?BR , , 1 

^ = BR(M^/.,)f;— , (22) 

where BR(M — > Ivi) is the total branching ratio of the corresponding non radiative 
decay. Finally, the total branching ratio BR is obtained by integrating Eq.(22) in the full 
kinematical range as follows 

BR^ [dx [dX (23) 
J J dx dX 

Q<x<l-ri, -^-^ < A < 1 . (24) 
1 — X 

In case in which kinematical cuts (x™°, and A™"^) should be apphed, the minima of 
integrations should be replaced as 

x"""" <x<l-ri, max jA'"^", ^-^| < A < 1 . (25) 



7 




Figure 2: Allowed helicity (in red) configurations of 7, and e"*" for 7r~^ — > e''~i^e7 decay in Tr"*" rest 
frame, figures (a), (b), (c), when all momenta (in blue) are aligned on the same axis. Direction 
of photon momentum is fixed by convention. Figure (d) corresponds to the non radiative decay 
vr'*' e^Uf.. Analogous spin configurations hold for the corresponding decays as well. 

Let us now consider the positron decay mode, where, as a good approximation, the 
lepton mass can be neglected in comparison to the pion one. A remarkable aspect of the 
results in Eq.(19), is that in the limit in which A — and x — 1, which corresponds 
to the emission of low energy positron and hard photons at relative small angles in the 
meson rest frame, the contribution proportional to /^(x. A) and to /^^j,(x. A) distributions 
dominates in the decay rate. In other words, hard photons will be mainly produced with 
left-handed polarizations. This behavior, as will be shown in more details in section 4, is 
just a consequence of the V-A nature of weak interactions and of the angular momentum 
conservation. This can be easily understood as follows. In the vr"*" rest frame, neglecting 
the lepton mass, we have x = 2E^/mM, and A = Ee/mM{l — cos^^e), where 6^e is the 
angle between positron and photon momenta. Let us consider the kinematical region in 
which A — s> 0, which corresponds to and/or 0. Due to the conservation 

of total momentum, and to the fact that \'Pe\/^M ^ 1 and 9^^ 0, the neutrino must 
be emitted in this region almost backward with respect to the photon direction, therefore 
final momenta are almost aligned on the same axis. This configuration is shown in Fig. 2a, 
where all momenta are chosen to be aligned on the same axis. One then could easily check 
the conservation of the spin {Sx) along the direction of the photon momentum which in 
Fig. 2 is set by convention on the negative X-axis. Since neutrino is a purely left-handed 
state, its spin projection along X-axis would be Sx{i^i) = +1/2. As a consequence of 
the angular momentum conservation {Sx{tt) = for pion), the photon must also be 
left-handed polarized giving Sxi'j) = — 1- In this case the positron, whose momentum 
is parallel to the one of the photon, must be right-handed in order to satisfy the total 
sum Sxi'j) + Sxi^e) + Sx{e~^) = 0. Notice that, in this particular kinematical limit 
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A ^ 0, photons with right-handed polarization would be suppressed, since the total sum 
of spins along X-axis would give in that case Sxil) + Sx{i^i) =3/2 thus spoiling angular 
momentum conservation. It is worth noticing that also in the case of 7r+ fj.^u^ decay 
mode, where the muon mass cannot be neglected in comparison to the one of the pion, 
the left-handed photon amplitude still dominates for high energy photons. This fact can 
be explained as follows. When the photon energy approaches its maximum value, being 
neutrino massless, in order to conserve total momentum, the production of the at rest 
it is favored. In this case the momentum of the neutrino should be opposite to one of the 
photon. As explained above, for this kinematical configuration, the photon is therefore 
favored to be produced as left-handed in order to conserve total angular momentum. 
Same considerations apply to the corresponding kaon decays as well. 

Another interesting case is the one in which the photon energy tends to zero, namely 
X ^ 0. In this singular kinematical region, one should expect soft photons to behave as 
scalar particles, carrying no spin. Then in this limit both the left-handed or right-handed 
distributions should tend to the same value, as indeed can be verified by performing the 
limit ,T ^ on the density distributions in Eq.(19). In the following we will show how 
this property could be relevant in order to define an observable which is free from infrared 
[E^ — > 0) singularity, namely the photon polarization asymmetry. 



3 The polarized radiative muon decay 

Here we analyze the radiative muon decay 

IJ,~{p) v^{qi) Pe(?2) e~(pe) l{k) 



(26) 



in which both photon and electron final states are polarized. The corresponding Feynman 
diagrams for this process are shown in Fig. 3, where p, p^, qi^2, k are the corresponding 
momenta. This decay is obtained from the non radiative one — > i/^PgC", by simply 
attaching the photon to the muon and electron external lines. Due to the V-A nature of 
weak interactions and a simple Fierz rearrangement, the square modulus of the polarized 
amplitude can be factorized as follows [4, 5, 15] 



,Ae)|2 



Gl 



(27) 



where M"(Ae, A-y) corresponds to the 0{q) amplitude in which photon is either radiated 
off the electron or off the muon, and A^"^ to the neutrino amplitude respectively 



M'^(A^,Ae) 
TV" 



+ ^ + me 

(Pe + ky 



-.it + ll 



m 



{Pn - kf 



mi > 



Uv{qi)llv^{q2) , 



u^{p) e*g{k, A^) 
(28) 
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(a) (b) 



Figure 3: Feynman diagrams for /i ^ z/^feC 7 decay. 

where u^{p), Ue{pe,K), Uu{qi,2) correspond to the muon, electron, and neutrinos four- 
spinors in momentum space respectively, with e and the corresponding helicities, 
and 'y2/R = (1 -F 75) 7"- Due to the factorization property of the amplitude in Eq.(27), 
one can easily calculate the sum over spins and the integral in phase space of neutrinos. 
At this purpose it is convenient to introduce the following tensor A^^^^ 

qi S" ( 
2 E\ 2 E2 



1 2 spins 

By making use of Lorentz covariance, one easily gets [15] 

A^"^ = Y ((P - Pe - kTiv -Pe- kY - g'-^ip -Pe- kf) • (30) 

where £^1^2 are the neutrinos energies. In order to describe the kinematic of muon radiative 
decay, we introduce the following independent variables 

2p-k 2p-pe 2k -pe 
x = — 5-, y = 5-, z = — , (31) 

mf, mi mi 

^ fj, fj, 

where m^ is the muon mass. In terms of these variables, the differential decay width, 
normalized to its tree-level non-radiative decay Fq, is given by 

dxdydz ~ ~2^ ^ Im^ ^' ^ 

Now we provide the expressions for the differential decay width in the rest frame of the 
muon, at fixed helicities of electron (Ag) and photon (A^), where as in previous section L 
and R symbols correspond to Ae,^ = — 1 and Xe^-y = 1 respectively. In particular, in the 
muon rest frame one has 

a;=^, ?/=— , z = ^{y-AcCos9) , (33) 
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where E^, arc the energies of electron and photon, cos^^ the angle between their 3- 
momenta, and Ag = \/y'^ — 4:r, with r = mg/m^. The allowed kinematical ranges for the 
above variables are 

< ^ < 2 f _ ^t-^"^ n ] . 2Vf<y<l+r, (34) 

\2 — y + cos J 

while — 1 < cos 6 < 1. Notice that the upper limit of x depends on cos 6. It is therefore not 
possible to perform the integration on cos^^ first. After the ^-integration the dependence 
on cos^^ it is also quite complicated. In the present analysis we are mainly interested in 
analyzing the structure of the leading logarithmic terms absorbing the regularized infrared 
and coUinear singularities as well as the one of the finite terms for the right-handed electron 
rate in the mg — > limit. For this purpose it is convenient to choose a particular region 
of the phase space where the analytical calculations are further simplified provided that, 
on the same time, all the leading logarithmic terms are preserved. As shown in Ref.[5], a 
suitable choice consists in taking the upper limit of x evaluated at cos 9 — 1, corresponding 
to its minimum value, as follows 

< ^ < 2 ( ^^^J'~^J , 2V^<y<l+r, (35) 

or equivalently 

2 \/f < y < ^-^^^^^^ , 0<x<l-y/^. (36) 

1 — X 

In this way, the integrals on x and cos 9 can be exchanged, giving a consistent simplification 
of the analytical integrations. 

For comparison, we will also provide the analytical expressions for the ^/-distributions 
and the total rates obtained by integrating over the full phase-space, but in the approx- 
imation of neglecting terms of 0{r). The corresponding analytical results at any order 
in r, will be presented elsewhere [16]. Regarding the corresponding numerical results, as 
shown in section 4, these are obtained by integrating the exact expression of the matrix 
density (given in Appendix B) at any order in r and over the full phase space. 

In the muon rest frame, the differential decay width is given by 

1 (ir^-^-^'-^^) air _ _ 1 

K {go + go) + K {gi + gi) , (37) 



Fo dx dy d cos 9 Stt x z'^ 

where the exact expressions at any order in r of the functions go,i and go,i, which depend 
on X, y, z, are provided in Appendix B. 

Notice that the r independent terms in the functions gfo,i and ^o,i ai'c proportional to 
z, partly compensating the 1/z'^ in front of the right hand side of Eq.(37). This is not 
true, however, for the r dependent terms which leave the distribution to be proportional 
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to Xj . These terms generate a singular behavior in the r ^ hmit for the distribution 
rate of the right-handed polarized electron, as it is in the analogous case of pion decay. 
Indeed, if the me — > limit is taken after integrating over cos ^, due to the property that 
/ d cos 9^2 oc 1/r, terms proportional to r/z'^ lead to a non- vanishing contribution in 
the integrated width. By taking into account the electron mass effects, the (polarized) 
integrated rate distributions (F^es'"'^'''') restricted range in Eq.(35) are given by 



1 dVi^j'^^^ ^ ^ 

Gq + Go + Ae {Gi + Gi) , (38) 



res ^ 



Fq dx dy 24:71 x 

where the expressions for the functions Go,i and (5o,i depending on x, y and r variables, 
are reported in Appendix B. 

The same phenomenon appearing in the meson decay for the right-handed electron 
[3] , it is also manifest here as a discontinuity in the electron mass. In particular, for right- 
handed electrons, the integrated rate in cos^ does not vanish in the nie — > limit, as 
one should expect from the massless theory. This discontinuity, firstly noticed in Ref.[6], 
can be associated to the axial anomaly [3, 8, 9] according to the interpretation of the axial 
anomaly given by Dolgov and Zakharov [7]. This anomalous behavior can be easily seen 
from the r — > limit of the functions Gq + Gi and Gq + Gi reported in Appendix B. For 
this purpose, we will provide below the expressions for the polarized differential decay 
width in the nie limit, in particular: 

lim — = ^—{-3(x-3)(x^-2) + n2 + x(9 + x(2x-5)))y 



+ 



(2x + 2y-3) (3 log(r)-6 1og(y))} 



hm r^— = 

'•-^o To dxdy 

lim-^^^%^ = ^-{4a;3(l + log(r)+|/)+2a; (|/-l)|/(l2 + 6 1og(r)+i/ 
^•^0 To dxdy 27i x ^ ' ^ 

+ 2 (2 + log(r)) (2 7/ - 3) + x^ (6 log(r) (2 y - 1) 
+ y (16 + 5y)-6) -4(x + y)' (2x + 2y-3) log(y)} 

lim — = - X {?>-2x-2y) , (39) 

where the log(r) terms are retained in order to regularize the coUinear divergences. 

The lepton is intrinsically left-handed, due to the nature of the coupling and parity 
violation. However, a final right-handed electron can also be produced with a sizeable 
rate in the limit r — > [5, 8]. As we can see from Eqs.(39), in the r — > hmit the photon is 
purely left-handed polarized when final electron is right-handed, as expected from angular 
momentum conservation. 
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QED obeys parity conservation and, therefore, the photon polarization is a "mea- 
sure" of parity violation. In the limit x ^ however the photon should behave as if its 
wavelength does not any longer resolve the process itself, leaving to a spin decoupling 
phenomenon as already observed in the radiative meson decay. Therefore, in this case, 
the soft photon does not take part to the angular momentum conservation of the whole 
process. Then, as in the meson case, in the a; — > limit the partial amplitudes corre- 
sponding to the distributions RL and LL, first and third above, tend to the same limit 
as well as the RR and LR ones. This property again shows the soft photon decoupling 
from any spin-related process. 

Now we provide the analytical results for the differential distribution in the electron 
energy, in the approximation of neglecting terms of C(r), obtained after integrating over x 
the distributions in Eq.(39). Since the total integral in x contains the well known infrared 
divergence when x — > 0, due to the emission of soft photons, we should provide integrated 
results by fixing a cut in the photon energy, namely Xq, corresponding to the experimental 
energy resolution of photon detector. Then, in the r — > limit, the kinematical range of 
xisxo<x<l — y. After integrating over the x range, and by retaining only the leading 
terms in xq and r, the result is 

^ = log(a;o)-log(l-|/) (3-2y)(2 + log(r)-2 1og(|/)) 

+ ^ (1 - y) (57 + 36 log(r) + 28 y + + 4 - 72 log(y)) } 

1 dr(5^) 



To dy 

1 (ir(£'^) 



dy 



= - |[log(xo)-log(l-y)ly2 (3 - 2 y) (2 + log(r) - 2 log(y)) 

- ^(y-lf (l0 + 96y + 5y2 + 2 1og(r) (5 + 22y) 

- 4 (5 + 22y) log(y))} 



To dy 



-(1 



yf (5 - 2y) . 



(40) 



Notice that the coefficient of the term proportional to log(xo) in Eqs.(40), should cancel 
the same term appearing in the one-loop corrections to the non-radiative Born decay, 
as shown in section 7. Here we would like to stress that the coefficients of the terms 
proportional to log(a;o), appearing only in the expressions of (i?, V) and (L, V) in Eqs.(40), 
are the same for both Left- and Right-handed photon contributions. This, again, shows 
the property that the photon spin must decouple in the infrared limit. 

In the zero lepton mass limit, the kinematical range of y are now < y < 1 — and 
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it is easy to check that the electron energy distribution vanishes at the end points. As 
a cross check of our results we integrate over y the non-vanishing expressions above and 
obtain 

r(f,'^) 23\ 17 , , , TT^ 10399. 

r(£.^) ar (I ,^5\ 17 17-. 

= - { log(-) (2 ^°S(^o) + 24 j + 12 - "6" + 18/ ' 



To 

Y{l,r) 



a 



To 47r 



(41) 



Finally, the total width for the radiative muon decay, integrated over the restricted phase 
space in Eq.(35), summed over all polarizations is 



afi /xA / N 7\ 17 , , , TT^ 13409 
- { log(r) (log(.„) + + - log(..) - y + ^} . (42) 



where terms of order 0{r) and 0{xq) were neglected ^. 

For comparison, we report below the ^-distributions integrated over the full phase 
space. As in Eqs.(40), we use the approximation of neglecting terms of order 0{r) and 
O{xo). In particular, for the polarized differential rates in the electron energy, we have 

1 rfr(^^'^«) 1 dPli-"^-) a 



+ ^A(^-^«), (43) 
To dy To ay tt 

where the additional terms arising from the extra phase-space integration, are 

given by 

A(^'^) = ^{-30t/+3y2(7 + 3L3(y))-t/^(37 + 6L3(2/)) + 27y^-3y^ + 4/} 
+ ^ log(l-|/) {-5 + 6|/ + 3|/' (-3 + 2|/) \og{y)-y^] , 

/^iL,L) ^ ^{32y-{95-6Ls)y' + {46-ALs{y))y' + 5y^} 

+ ^{\og{l-y) (5-24y + 30y^-lly^-3y^ {3 - 2 y) \og{y)) } , 
A(^'«) = (44) 



^Here we would like to stress that this expression agrees with the corresponding one reported in 
Ref.[5], but differs from the old results in Refs.[17, 18]. We remind here that the total width calculated 
in Refs.[17, 18] is obtained by integrating over the full phase space. Therefore, the coefficients of the 
logarithmic terms coincide with the corresponding ones in Refs.[17, 18], as expected since both infrared 
and coUinear singularities arc included in the phase space region of Eq.(35), Therefore, the total width in 
Refs.[17, 18] will differ with respect to Eq.(42) by finite non logarithmic terms in the xo ^ and r — » 
limits. 
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with L^i^y) = TT^ — 6Li2(l — y). As shown in Eqs.(44), only the LL and RL distributions 
get an extra contribution which is non- vanishing in the r ^ limit, while for the cor- 
responding RL and RR ones this is of order 0{r). This is due to the fact that in the 
radiative muon decay the right-handed electron is mainly produced at ^ ~ 0. Hence, 
regarding the right-handed-electron production, the maximum of the x range integration 
can be well approximated by a;™^^(cos^) ~ x'^'^^{cos9 — 1). This approximation, adopted 
in Ref.[5], corresponds to consider the restricted phase space in Eq.(35). 

Finally, by integrating the distributions in Eqs.(44) over the full range < y < 1, we 
obtain 



r(^'-^) a r , , , /I , , 23 
-TT = ;^l^°S^^)l2^°g^^°^ + 24 

- I log(r) log(xo) + 



+ 



To TT °' ^2 °' 24 

Then, the total rate summed over all polarizations is given by 



17 


log(a;o) - 




997 


12 


12 ' 


^ 288 


17 


log(xo) - 




133 






12 


12 ' 


^ 288 



(45) 



= - { log(r) (log{xo) + g j + y log(^o) - y + ^} • (46) 



r, 



which is in agreement with the previous result obtained in [17-19]. We would like to 
stress here that the structure of the leading logarithmic terms is also preserved in the 
polarized rates when the restricted phase-space-integration is considered, as can be seen 
by comparing Eqs.(41) and (45). 



4 Distributions and polarization asymmetries 

In this section we present the numerical results for the distributions of branching ratios in 
the photon and electron energies. In both cases we will sum over the fermion polarizations, 
leaving fixed only the photon polarizations. For this purpose, it is very useful to introduce 
also an observable which provides a direct measurement of the amount of parity violation 
in the weak decays, namely the distribution of photon polarization asymmetry A-j,, defined 
as follows 

dA, ^ d^(BRL)-dg(BRR) 

dC d5(BRL) + d4(BRR) ^ ^ 

where d^(BRL,R) = '^^^'"'"^ stands for the differential branching ratio (BR) in ^ = {x,y}, 
where x = 2Ej/M and y = 2Ei/M in the rest frame of the decaying particle of mass M. 
Labels L and R indicate left- and right-handed photon polarizations respectively. Here we 
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would like to stress that —rr is a finite quantity, free from infrared divergences. Indeed, 
when the photon energy goes to zero, the distribution tends to zero, since 

lini {pi^{x, y) - pr(x, y)} 0{x) and lini {pi^{x, y) + pr(x, y)} log(x) ,(48) 

where piji{x,y) indicates a generic Dalitz plot distribution for the polarized decay with 
left- (L) or right-handed (R) photons. Therefore, the total integral of Eq.(47) is a finite 
and universal quantity. It does not depend on the photon energy resolution of the detector, 
and provides a direct measure of parity violation. Moreover, being quite sensitive to 
the hadronic structure of radiative meson decays, it is also an useful tool for accurate 
measurements of V and A form factors. 

In the following sections 4.1, 4.2, and 4.3, we will show our results separately for the 
case of pion, kaon, and muon decays respectively. Let us start with pion decay. 



4.1 Radiative 7r+ decays 

Here we report the numerical results obtained for the distributions of BRs and asymme- 
tries for the case of radiative pion decays tt"^ — > e'^Ue^ and tt"*" — > As shown in 
section 2, the corresponding amplitudes contain only two free parameters which enter in 
the hadronic structure dependent terms (SD), that is V and A form factors. However, 
being V and A non-perturbative hadronic quantities, they cannot be evaluated in QCD 
perturbation theory. An alternative approach like the one of effective field theories as, for 
instance, Chiral Perturbation Theory (ChPT), or the lattice QCD, should be employed. 
On the other hand, V and A could be directly measured by experiments [10,20-22]. 

These form factors are not constant over the allowed phase space. Nevertheless, in 
radiative pion decays, the momentum dependence in and A{W^), parametrized by 

W'^ = (l — x) m^, is expected to be very small, not exceeding a few per cent of the allowed 
phase space. This expectation is also supported by ChPT, since at the leading order in 
ChPT V and A are constant. Recent calculations at next-to-leading order in ChPT [23], 
which included terms up to 0{p^), where p indicates a generic momentum involved in the 
decay, show a mild dependence on momenta, confirming the above expectations. In our 
analysis, we will assume V and A to be constant in the full kinematical range^. Now we 
summarize the present status of form factors determination in pion decay. 

The vectorial form factor V can be extracted in a model independent way from 7r° 
77. By using the conservation of vectorial current (CVC) hypothesis, one can relate the 

^Taking into account the effect of momentum dependence in the form factors goes beyond the purpose 
of the present work, since we are mainly interested in analyzing the dependence of photon polarization 

asymmetries by the photon and electron energies. For more accurate predictions of BRs and asymmetries, 
these effects should be included, especially in the kaon decay where they are expected to be sizeable. 
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vectorial form factor to the lifetime of the neutral pion [24] 



1 /2r(7rO ^ 77) 
V\ = -W — ^ = 0.0259 ± 0.0005 , 



(49) 



where r(7r° — > 77) is the total width of 7r° — > 77 decay and V is assumed constant. On the 
other hand, the axial A form factor can be measured via the ratio 7 = V/A. In previous 
experiments [20], using the stopped pion technique, the radiative pion decay has been 
measured in a limited phase space region where V + A contributions dominate, leaving 
to an ambiguity on the sign of 7. In more recent experiments [21,22], the investigated 
larger portion of the phase space allowed to determine the sign of 7 as well, which has 
also been confirmed by the tt"'" — >■ e"'"i/e^e~ measurement [25]. 

The most recent measurements of radiative pion decay, using the stopped pion tech- 
nique, has been performed by the PIBETA collaboration with a good accuracy [10] . There, 
the CVC hypothesis has been used for the determination of 7. The preliminary results 
of PIBETA experiment indicate a deficit of events in the observed tt — > ei/7 decay [10], 
suggesting for a new tcnsorial four-fermion interaction beyond the V-A theory. An analo- 
gous effect was first observed in a previous experiment at ISTRA facility in early 90s [22] , 
in which pion decays where studied in flight. We will return on this point in section 6, 
where the potential role of new tensorial couplings, suggested in order to accommodate 
experimental data, will be discussed. 

The results contained in this section have been obtained by using for 7 the central 
value of the best CVC fit reported by the PIBETA experiment, namely 



which is also consistent with predictions in ChPT. 

In Fig. 4 we show the distributions of BR for pion decay in electron channel. In par- 
ticular, dBR/dx and dBR/dy are reported in the "top" and "bottom" plots respectively^. 
In the plots we have integrated the phase space over A and x respectively. Results are 
obtained by imposing kinematial cuts on A or x, as indicated in the plots. Looking at the 
dBK/dx distributions in Fig. 4, the general behavior emerging from these results is the 
following. When cuts on A are relaxed, the left-handed photon polarizations dominate in 
all range of values oi x > 0.2, corresponding to photon energies Ej > 14 MeV. On the 
other hand, the contribution of right-handed photons can be increased by imposing larger 
cuts on A, as can be seen by comparing left-top and right-top plots in Fig.4. For example, 
by requiring that A > 0.3, right-handed polarizations could dominate in the region of 
hard photons 0.6 < x < 0.9. These results can be explained by using angular momen- 
tum conservation. When cuts on A are relaxed, the main contribution to the integral in 

^In the following, in all the plots involving the distributions of BRs, green and red curves correspond 
to left- (L) and right-handed (R) photon polarizations respectively, as also indicated in each curve. 
Unpolarized decays are drawn as blue curves, with L-l-R label associated to them. 



7 = 0.443 ±0.015, 



with 



V = 0.0259 



(50) 
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Fi gure 4: The photon energy spectrum — ^ versus x (top plots) and electron energy spectrum 
'^^y^ versus y (bottom plots), for pion decay 7r+ — > 1^66+7. The labels L and R attached 
to the curves indicate pure left-handed and right-handed photon polarizations contributions 
respectively, while L + R correspond to the sum. Kinematical cuts A > 0.1 (top-left), A > 0.3 
(top-right) and x > 0.3 (bottom- left), x > 0.6 (bottom-right) are applied respectively. 
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dX comes from the region of low A, where the IB effects dominate with respect to SD 
terms. Low values of A should correspond to small angles between photons and e+, but 
could also correspond to low positron energies. In the former case, neutrinos are likely 
to be produced with opposite direction with respect to the photon momentum, in order 
to compensate for the missing momentum in the pion rest frame. Since neutrinos are 
always left-handed polarized, photons must be left-handed as well, as required by angular 
momentum conservation. The spin configuration for this case is shown in Fig. 2a. On the 
other hand, small values of A could also correspond, in the latter case, to the spin con- 
figuration shown in Fig. 2c, where positron and photon are backward. There, the photon 
should be mainly emitted from the tt"*" line, leaving positron and photon both left-handed 
polarized. As we will show in section 5, after integrating over x with cuts x > 0.3, the 
dominant effect will be given by this last configuration. 

On the contrary, when cuts on A are very large, IB effects are reduced and SD terms 
become sizeable. In this case, the positron is mainly produced right-handed, due to the 
fact that SD terms are not chiral suppressed, neutrino momentum is favored to be directed 
forward with respect to the photon one, leading to a right-handed photon as shown in 
Fig. 2b. However, as we can see from the top plots in Fig. 4, there is a region of large x 
where left-handed photon contributions are also sizeable, in particular for 0.9 < x < 1. 
This peculiar behavior in the end point region of photon energy can be explained as 
follows. When photon energy approaches its maximum, positrons start to be produced 
almost at rest, if A is small. Then, in order to conserve the total momentum, neutrino 
should be mainly emitted backward with respect to the photon direction, see Figs. 2a 
and 2c, leading to left-handed photon polarizations as required by angular momentum 
conservation. However, we should stress that, depending on the cuts on A, the right- 
handed photon polarization could dominate even near the end-point region of photon 
energy.^ 

In the bottom plots of Fig.4, we report the BRs distributions on positron energy versus 
y, for two representative choices of cuts, namely x > 0.3 (left-plot) and x > 0.6 (right- 
plot). As we can see from these results, the left-handed photon polarizations dominate 
in the region y < 0.5, while the gap between left-handed and right-handed contributions 
increases by using stronger cuts on the photon energies. This behavior can be roughly 
understood as follows. At fixed positron energy, the larger the photon energy is the 
more the neutrinos are produced parallel and backward to the photon direction, in order 
to conserve total momentum. Therefore, as explained above, conservation of angular 
momentum favors in this case the left-handed photon polarizations. However, at the end 
point of positron energy, when cuts on a; > 0.3 and x > 0.6 are imposed, the scenario 
could be reversed. As a consequence of the strong cuts on x, the IB contribution can be 

^Even if it is not shown in the plot, at the end point x = 1 — re the L curves (as well as R and L+R 
ones ) corresponding to the dBR/dx distributions in Fig.4 vanish. However, the L curves go to zero more 
slowly than the corresponding R ones. 
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Figure 5: As in Fig. 4, but for pion decay vr"*" f^//+7, and with kinematical cuts A > 0.6 
(top-left), A > 0.8 (top-right) and a; > 0.1 (bottom-left), x > 0.2 (bottom-right). 



made very small, and near the region oiy = Umax-, the SD terms should dominate favoring 
the production of a right-handed positron. Clearly, when photon and positron are both 
very energetic they tend to be emitted with opposite direction in order to conserve total 
momentum, approaching, in the case of a right-handed positron, to the spin configuration 
in Fig. 2b. Therefore, due to angular momentum conservation, photons are mainly right- 
handed in the region y > 0.8 and x > 0.3. Here we would like to stress that the relative 
gap between left- and right-handed contributions of hard photons, near the region y > 0.8, 
should be very sensitive to the values of hadronic form factors. 

In Fig. 5 we show the corresponding results for the pion decay in the muon channel. 
In this case we can see that the left-handed photon polarizations always dominate over 
the entire phase space, while right-handed ones are quite suppressed. Notice that, being 
the muon mass very close to the pion one, the IB contribution is not chiral suppressed as 
in the electron channel and it is larger than the SD one almost over all the allowed phase 
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space. This implies that /i+ is mainly produced with left-handed polarization. Moreover, 
due to the fact that here the minimum allowed value of A (for a; > 2) is Amin ~ 0.7, 
the and photons are mainly produced at large angles. Then, if left-handed /i'^ tends 
to be produced backward with respect to photon momentum, this last one must be also 
left-handed in order to conserve total angular momentum, as shown in Fig. 2c. 

Finally, in Fig.4, we present our results for the x- and ^-distributions of the photon 
polarization asymmetry, as defined in Eq.(47). In particular, in the top and bottom plots 
we report the results for the dA^/dx and dA^/dy respectively for several kinematical 
cuts, while the left and right plots correspond to the electron and muon channel decays 
respectively. A general property of these results is that the x- and ^-distributions of 
asymmetry vanish at a; = and y = 1 + ri respectively. This is a consequence of the fact 
that when the photon energy is approaching to zero, the polarized photon densities of 
Dalitz plot tend to the same limit, due to the spin-decoupling property of soft photons, as 
discussed in section 2. A remarkable aspect of these results is that, in the electron decay 
channel, the dA^/dx distribution becomes negative for some particular choices of cuts. 
Analogously, the same effect can be achieved on the y-distribution by increasing cuts on 
the photon energy. On the contrary, in the muon channel, the corresponding asymmetry 
is always positive, as can be seen in the plots to the right in Fig.6. In conclusion, we 
would like to emphasize that the position of the zeros of photon polarization asymmetry 
is particularly sensitive to the effects of the SD terms. This property could suggest a new 
experimental way for obtaining more precise measurements of form factors. 

4.2 Radiative decays 

In analogy with the radiative pion decays, we analyze here the corresponding ones in the 
kaon sector, in particular —>■ e^v^. 7 and — > \i^v^ 7. The expressions of amplitudes 
in terms of decay constants, masses and form factors remain formally the same as in the 
pion decay. However, the kaon electromagnetic form factors, as well as the decay constant 
jxi and the ratios Tg, r^^ between leptons and kaon mass, are quite different from the pion 
case. As we will see in the following, these differences will sizeably affect the shape of 
distributions and asymmetries with respect to the corresponding pion decay. 

The most recent measurements of V and A form factors have been performed by 
the E787 collaboration [26] through radiative decay fi'^u^'-f. In particular, the 

absolute value of V+A has been determined finding |^+^| = 0.165±0.007±0.011, while a 
limit on —0.24 <V — A< 0.04 has been set at 90% C.L. These results have been obtained 
by assuming constant form factors. The |V^+^| measurements are consistent with previous 
results on X — > e'^Ug^, but they disagree by almost 2 standard deviations with respect 
to predictions from leading order in ChPT [12]. We recall here that the evaluation of 
form factors starts at one loop in ChPT expansion, that is at C(p^). At this order 
the chiral prediction, as for the pion case, gives constant form factors. The momentum 



21 





A 
I 

> 

A 
I 

+ 
< 

■o 



0.8 



0.6 



0.4 



0.2 



' \ ' I 






0.2 




0.1 




X > 



1.5 1.51 1.52 1.53 1.54 1.55 1.56 1.57 



Figure 6: The differential asymmetry versus x (top) and versus y (bottom), witli 
kinematical cuts A > 0, 0.1, 0.2, 0.3 (top-left), A > 0, 0.8 (top-right) and y > 0, 0.2, 0.4, 0.6 
(bottom-left), y > 0, 0.1, 0.2, 0.3 (bottom-right) respectively. 
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dependence of the form factors starts then at the next-to- leading order, that is 0{p^), and 
it is expected to be larger than in pion case, due to sizeable effects of resonances exchange 
[27]. In particular, by considering only a particular class of diagrams where vector and 
axial-vector resonances are exchanged, the form factors can be parametrized as 

1 - W^/mj^* 1 - W^/mj^^ 

where — m'j^{l — x), and the masses m|-* and mj^_^ correspond to vector and axial- 
vector resonances. Then, in order to minimize the effects of resonance exchange, large 
x-regions should be considered since — > when x — > 1, while low x-regions may 
be used to explore the dependence of V and A. The 0(p^) contributions, based on 
SU{3) X SU{3) symmetry in ChPT, has been recently calculated in Ref.[28]. Significant 
deviations of order of 10-20 % have been found on V and A with respect to the leading 
order calculation. Moreover, while the vectorial form factor is quite sensitive to photon 
energies, the axial one shows only a modest effect [28] . 

As for the pion decay, in order to simplify our analysis, we will not take into account 
the momentum dependence in V and A. Then, consistently, we will take the V, A pre- 
dictions at the leading order in ChPT, re-absorbing the missing NLO contributions in 
the theoretical uncertainty. In particular, our results are obtained by using the following 
values [12] 

V + A^ -0.137, y- A = -0.052. (52) 

In Fig. 7 we show the x— (top) and y— (bottom) distributions for the — >• e'^u^j 
decay. The general trend emerging from these results is that in kaon decay, contrary to 
the pion case, the right-handed photon production dominates over the lef-handed one, 
already for moderate cuts A > 0.1, as can be checked by comparing results between Figs. 4 
and 7 with the same cuts on x and A. This result can be roughly understood as follows. 
The IB contribution in the radiative decay is more "chiral" suppressed with respect 
to the corresponding tt"*" due to the fact that rriK — 4:771-,^. Then, for the same values of 
X and A, the IB effects in pion decay will be larger than in the corresponding kaon one. 
For instance, while the IB contributions in pion decay are still sizeable after cuts A > 0.1 
and X > 2 have been imposed, in if + decay these same cuts dramatically reduce the IB 
effects in favor of SD contributions. As already explained in section 4.1, when the photon 
is produced from the SD terms it is mainly right-handed polarized. In conclusion, the 
dBR/dx distributions in the top-plots of Fig.7, for x > 0.1, shows the same behavior of 
the corresponding one in pion decay in the region 0.6 < x < 0.9 and A > 0.3, where 
the right-handed photon contributions are enhanced. Moreover, as we can see from the 
top-plots in Fig.7, these curves have a maximum (for 0.1 < a; < 1 and 0.1 < A < 0.3) 
around x ~ 0.75. In the bottom-plots of Fig.7 we report the analogous results for the 
positron energy distributions. As we can see, when the y > 0.5 and a; > 0.1 cuts are 
imposed, the right-handed photon gives the dominant effect. 
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In Fig. 8 results for BR distributions in x and y are shown for the — > decay. 
The hierarchy between L and R curves, and their shapes, are similar to the corresponding 
ones of 7r+ /W^t'/iT ( see Fig. 7). Notice that the available ranges of x and y are larger for 
the kaon decay in the muon channel with respect to the corresponding pion decay, due to 
more available phase space of the former. The same considerations regarding the shapes 
of X- and ^-distributions of pion decay should hold here as well. As can be seen from 
these results, also in —>■ /i"^z/^7 decay the left-handed photon polarization gives the 
dominant effect for x > 0.2 and A > 0.3 or analogously for x > 0.3 and y > 0.4, as shown 
in the bottom plots for the y distribution. Numerical results for the total contribution 
L-l-R, are consistent with the corresponding ones in Ref. [11]. 

In Fig. 9 we show our results for the distribution of asymmetries in the kaon decays. As 
we can see from left-top plots in Fig.9, the dA^/dx distribution for K'^ e^Vel is always 
negative in all range x > 0.1, and tends to a constant value for x > 0.3 — 0.4, depending 
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(a) (b) (c) (d) (e) 

Figure 10: Some helicity (in red) configurations of 7, and e~ for [i~ — > v^Vf-e^^ decay in 

\x~ rest frame, figures (c)-(e), wlien all momenta (in blue) are aligned on the same axis. Figures 
(a) and (b) correspond to the non-radiative decay [i~ — > VyVe- Direction of electron momentum 
in (a)-(b), as well as photon momentum in (c)-(e) diagrams, is fixed by convention. 

on the applied cuts on A. In particular, when A > 0.6, the dA^/dx already approaches 
its minimum value for x > 0.2. By relaxing the constraints on A, we can see that the 
dA^/dx distribution could have a zero at a; ~ 0.2, and analogously dA^/dy at ?/ ^ 0.5. 
On the right-plots we present the corresponding results for the I^^^plI decay, and 

for some representative choices of cuts. As we can see, the dA^/dx is more sensitive to 
cuts on A than the corresponding one in vr"*" — >■ Analogous considerations hold 

for the dA^/dy distribution as well. In conclusion, the photon polarization asymmetry 
for radiative meson decays in muon channel, is always positive, vanishing only at the end 
point X = or analogously y = I, due to the spin decouphng property of the soft photon. 

4.3 Radiative /i decay 

In this section we will discuss the numerical results for the radiative muon decay fi^ —>■ 
v^j^e'Vel- As shown in section 3, this process is described by the leptonic Fermi interaction, 
where the photon is attached to external legs of muon and electron, see Fig. 3. Being a 
pure leptonic process, its decay rate can be calculated with high accuracy in perturbation 
theory. In particular, the 1-loop QED corrections have been evaluated in Ref.[15] for the 
inclusive radiative muon decay, which corresponds to an accuracy of order 0{a'^) in the 
branching ratio. However, studies of polarized radiative muon decays have been recently 
published [4,5]. Also 1-loop radiative corrections have been included in the evaluation of 
the decay rate [4]. However, in these studies only the polarization of fermions has been 
considered. 

Our results for the polarized photon distributions of branching ratios are shown in 
Fig. 11. These results, as well as the corresponding ones in Figs. 12 and 15, have been 
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Figure 11: The photon energy {E^) spectrum '^^^^ versus x = 2E^/m^ (top) and electron 
energy (E-y) spectrum — versus y = 2Ee/m^ (bottom) for muon decay i>^U(,e^ with 

left-handed (L) and right-handed (R) photon polarizations, and for kinematical cuts y > 0.2 
(top-left), y > 0.6 (top-right) and x > 0.2 (bottom- left), x > 0.6 (bottom-right) respectively. 

obtained by integrating over the full the phase space and by taking into account the full 
r dependence. From Fig. 11 we can see that the main contribution to the radiative decay 
is provided by the left-handed photon polarization, while the right-handed one is quite 
suppressed and decreases by increasing the photon energy. This behavior, again, can be 
explained by using angular momentum conservation and parity violation. Notice that, 
due to the V-A nature of weak interactions, the electron is mainly produced left-handed 
polarized in the muon decay, and chirahty flip effects, needed to produce a right-handed 
electron, are always sub-leading, being proportional to the electron mass. Moreover, due 
to the fact that we are integrating over the final phase space of neutrinos, the analysis 
is strongly simplified. Indeed, after integration, the effect of neutrinos is re-absorbed 
in the tensor A^"^ appearing in Eqs.(29), (30). Notice that A^"^ is just a projector for 
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the four-momentum = —{Pe + k + p)" carried by the neutrinos pair, and it can be 
seen as the sum over polarization states of a massive particle of spin 1. In other words, 
regarding the spin content, the neutrinos pair behave as a spin-1 particle of mass Q^, 
having three polarization states. In the case of non-radiative muon decay, the allowed 
spin configurations in the muon rest frame are shown in Fig.lOa,b, where all momenta are 
aligned on the same axis X and by convention the electron momentum is chosen along 
the negative direction. As we can see, if the electron is left-handed (Jx = 1/2), the spin 
projection of neutrino anti-neutrino pair along the direction of their total momentum can 
be Jx = —1 or Jx = 0, but not Jx = +1, being the muon a spin 1/2 particle. 

Let us now consider the radiative decay, with photon emission from the electron line. 
It is known that, when hard photons are emitted parallel and forward to the electron 
momentum, they can flip the electron helicity, without paying any chiral mass suppression 
[3-8]. The helicity-fiip mechanism is illustrated below 

\ ^ \^ 





(a) (b) 

for incoming left-handed (a) and right-handed ept (b) electron by collinear photon 
bremsstrahlung. All momenta, indicated by blue arrows, are aligned on the same axis. 
Red arrows stand for the corresponding helicities and an electron mass insertion is under- 
stood. We remind here that the chiral suppression of the term ml appearing in the square 
modulus of the numerator due to the chirality flip, is compensated by the singular be- 
havior in ~ 1/ rrig appearing in the square modulus of the propagator for collinear photon 
emission. The corresponding spin configurations in this case are shown in Fig.lOd-e for 
the case of aligned momenta. However, as for the meson case, the contributions coming 
from the helicity-fiip transitions are always smaller with respect to the helicity conserving 
ones. The largest contributions should come from the photons emitted by the muon line. 
In this case the neutrinos and electron momenta are favored to be busted backward with 
respect to the photon momentum, in order to compensate for the missing momentum. 
The corresponding spin configuration, in the particular limit in which all momenta are 
aligned on the same axis, is shown in Fig. 10c. Since the favored spin of the u + u system is 
in this case Jx = 0, the photon must be necessarily left-handed polarized. This peculiar 
configuration should explain why the left-handed photon contribution is always dominant 
with respect to the right-handed one, leading to an increasing relative gap as the pho- 
ton energy increases. This seems to be the case since, as the photon energy approaches 
the soft region x ^ 0, the gap should decrease due to the spin decoupling property of 
soft photons. Results concerning the branching ratio distributions for the production of 
right-handed electron are shown in the next section. 
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Figure 12: As in Fig.ll, but for the asymmetry versus x (left) and versus y (right) 
and for kinematical cuts y and x > {0, 0.2 , 0.4 , 0.6} for left and right plots respectively. 

In Fig. 12 we show in the right (left) plots the dA^/dy (dA^/dx) distributions asymme- 
try for the radiative muon decay. We present our results for some representative choices 
of cuts, in particular y > {0, 0.2, 0.4, 0.6} and x > {0, 0.2, 0.4, 0.6} for the dA^/dx and 
dA^/dy respectively. As we can see from these results, the asymmetry in the muon case 
is always positive, as a consequence of the dominant left-handed photon contribution as 
discussed above. The shapes of asymmetries are quite similar to the corresponding ones 
in 7r+ — i> fi~^i'fj,'~f and — >■ /i"''i/^7, where the IB effects are larger than the SD terms. 
In particular, here the dA^/dx is monotonically increasing with x, while analogously the 
dA^/dy is monotonically decreasing. 

5 Energy spectra of the polarized positron/electron 

In this section we discuss the results for the distributions of BRs in the positron energy 
for the decays 7r+ e^Vel and K'^ —>■ e^Veli and analogously in the electron energy for 
the radiative muon decay, for both lepton and photons polarizations. As seen before, due 
to angular momentum conservation, the IB contributions favors a left-handed positron in 
the soft photon energy region. However, in the case of an hard photon produced collinear 
with e^, the IB favors a right-handed positron. This last effect is due to the helicity-flip 
mechanism discussed above. However, if one imposes larger cuts on x in order to reduce 
the IB effect, then the relative contribution of the SD term grows up. In particular, when 
the positron is totally generated from the SD terms, it is mainly right-handed polarized 
(see discussions in section 4.1), as can be seen from Eq.(75) in Appendix A, since the 
photon is emitted from the hadronic vertex. From these results one can easily check 
that in the region of ?/ ^ 1, the LR term in Eq.(75) tends to zero and survives only 
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Figure 13: The electron energy spectrum versus y for vr"*" — > fee"'"7, for left-handed (left 
plot) and right-handed electron polarizations (right plot), with photon energy cut x > 0.3. As 
in previous figures, the labels L and R labels attached to the curves indicate pure left-handed 
and right-handed photon polarizations contributions respectively, while L + R correspond to the 
sum. 



the RR one, corresponding to the production of both positron and photon right-handed 
polarized, as required by angular momentum conservation. In particular, for a generic 
meson M = if, vr, we have 

, ^ ^ ^iy±^< (1 _ + . (53) 

y^i ay zvr 48 n ^ ^ 

where x^ut is the cut onx, BR(7r+ e+z/) = 1.23x10"^ and BR(if + e+z/) = 1.55x10"^ 
[29]. Notice that for a right-handed electron, the contribution from the other polarizations 
vanishes in the massless lepton limit r/ — > and for y — > 1. 

In Fig. 13 we show the plots corresponding to the dSK/dy of pion decay, where on the 
left and right plots we report the case of left-handed (ej) and right-handed (e^) positron 
polarizations respectively. These results are obtained for x > 0.3. As we can see, the 
contribution of is always dominant with respect to one. This is because the IB effect 
is still large for x > 0.3, and so a left-handed positron is favored. Moreover, a peculiar 
aspect of these results is that for the production, the left-handed photon polarization 
dominates in all the positron energy range. On the other hand, in distribution, the 
right-handed photon gives the leading effect for y > 0.7, being totally induced by the SD 
terms. The maximum of dBK/dy for production, achieved at ?/ ^ 1, can be easily 
checked by using the approximated expression in Eq.(53). 

Analogous results for the kaon decay are shown in Fig. 14. Here, the situation is 
reversed with respect to the pion case. The contribution gives the leading effect in 
the total BR already for x > 0.3 and the right-handed photon polarization dominates 
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Figure 14: As in Fig.13, but for —>■ Ve^^l- 

for y > 0.5. On the other hand, in the distributions (see left plot in Fig. 14), the 
left-handed photon contribution provides the dominant effect. As already mentioned in 
section 4.2, these differences with pion case are mainly a consequence of the fact that the 
IB amplitude is more chiral suppressed in K —>■ eu^'y than in vr — > euej. 

Finally, in Fig. 15 we present the electron energy distributions for the muon case. In 
this case we imposed a cut x > 0.2 on the photon energy. We see that in both cr and ei 
distributions, the left-handed photon contribution always provides the dominant effect, 
being this configuration favored by angular momentum conservation. Moreover, in the 
case, the right-handed photon contribution is very tiny. As we explained in the previous 
section, these results are a consequence of the fact that in the radiative muon decay the 
electron is naturally produced left-handed due to the V-A theory. On the other hand, the 
contribution of e^j is mainly generated from the helicity flip mechanism of hard photons 
emitted from the electron line and it is a sub- leading effect. 

6 Tensorial couplings 

Here we analyze the dependence of the photon polarization asymmetry, in the radiative 
pion and kaon decays, induced by tensorial couplings. The aim of this study is motivated 
by the recent measurements of the radiative pion decay vr"*" [10], where a sig- 

niflcant discrepancy in the branching ratio, with respect to the SM predictions [10], has 
been observed. This anomaly might be interpreted as the effect of a centi-weak tensorial 
interaction beyond the V-A theory [14]. 

An analogous discrepancy was noticed long time ago at the ISTRA experiment [22], 
where radiative pions decays were studied in flight. In that experiment, the 7r+ — u^e^'y 
was investigated over a large phase space region, in particular 0.3 < x < 1 and 0.2 < A < 1. 
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Figure 15: The electron energy spectrum versus y for /i+ P^fee+, for left-handed (left 
plot) and right-handed electron polarizations (right plot), with photon energy cut x > 0.2. The 
curve in the right plot corresponding to right-handed photon (R), has been multiplied by a 
rescaling factor of 10^. 



The measured branching ratio = (1.6 ± 0.23) x 10~^ [22] was found significantly 
smaller than the expected one B^^ = (2.41 ± 0.07) x 10^'', based on the CVC hypothesis 
and V-A theory of SM. The fact that the measured number of events is less than expected, 
cannot be explained by a missing unknown background. This result was interpreted [14] 
as a possible indication of a tensorial quark-lepton interaction with coupling of order 10~^ 
in unity of Gp. In particular, the suggested new contribution to the effective Hamiltonian 
for AS = transitions is [14] 



^efr° = ^^0^ Vud [u(T,,{l - 75)rf] [ea,,{l - 75)^^e] + h.c. (54) 



where a^^, = 1/2 [7^, •y^] and fx a dimensionless coupling. Notice that tensorial interactions 
are not subject to the strong constraints coming from the non radiative decay tt — *• i^e 
(as, for instance, for the scalar interactions) simply because, the Lorentz covariance forces 
the hadronic matrix element (0| [ua^u{^ — 75)*^] l^r) to vanish. On the other hand, H^ff^^ 
can contribute to the amplitude (M-r) of the radiative decay 7r+ — >• I'eG^I as 

Mt = Kd Fre^* q'' [ea^,{l - 75)//^] • (55) 

The constant Ft can be related to fx in Eq.(54) by using low energy theorems and PCAC 
hypothesis [14], [30] 

F'T = \^fT (56) 
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where (/x) = (0|gg|0) is the vacuum expectation for the quark condensate and x is defined 
by [31] 



{0\qa^,q\e{k)) = e,x {0\qq\0) F^, , (57) 

with Cq the quark (g) electric charge, F^ty = i{k^ty — k^e^), and e^{k) the photon polariza- 
tion vector of momentum k. Then, the destructive interference between the SM and tenso- 
rial amphtudes accounts for the correct number of "missing" events observed at ISTRA if 
Ft = (5.6±1.7) X 10~^ [14], corresponding to atensorial coupling fr ~ (1.4±0.4) x 10~^ 
[30], where (//) = -(0.24 GeV)^ and x = -(5.7±0.6) GeV~^ values have been used [31]. 
This result is consistent with the limit < 0.095 (at 68% confidence level) coming from 
beta decay [14]. 

Clearly, if confirmed, this effect would be a clear signal of new physics. Indeed, in the 
SM, the tensorial coupling fj,^ is very small, being generated at two-loop level and chiral 
suppressed. In Ref . [30] the supersymmetric (SUSY) origin of /t has been analyzed. The 
leading SUSY contribution to /r, given by charginos and squarks exchanges in penguin 
and box diagrams, can be larger than SM one since it is induced at one-loop. However, 
present bounds on SUSY particle spectra do not allow fr to be larger than fx ^ 10~^, too 
small for the required value suggested in [14]. Moreover, there have been also criticisms 
about the consistency of such large tensorial couplings. In Ref. [32] it was pointed out 
that, due to QED corrections, an Jt of order of (9(10^^) might run in troubles. Indeed, 
the operator in Eq.(54) can mix under QED radiative corrections with a scalar operator, 
whose contribution is strongly constrained by n'^ e+z/g [32]. In particular, an upper 
bound on fx < 10~^ can be set by imposing the strong constraints on scalar interactions 
coming from tt'^ — > e'^Ue, which is two order of magnitude smaller than the required one 
[14]. However, more accurate analyses showed that it is possible to relax or even avoid the 
upper bound claimed in [32]. For instance, the simultaneous (fine-tuned) contributions 
of both tensorial and scalar interactions, as suggested by lepto-quark models [33], might 
relax the upper bound in [32] and thus the interpretation given in [14] cannot be regarded 
yet as ruled out. There is also an alternative solution, proposed in Ref. [34], where a 
modified tensorial interaction can formally avoid the mixing with scalar operator, while 
solving the ISTRA discrepancy. 

In Ref. [11], it was pointed out that an analogous effect might show up in the kaon 
sector. In particular, if the origin of fr is flavour independent, then a tensorial interaction 
of the same order is also expected in the AS — 1 transitions, leaving to a large anomaly 
in radiative kaon decays, easily detected at present and future kaon factories [11]. 

Recently, the PIBETA collaboration at Paul Scherrer Institute facility, has performed 
an accurate analysis of the tt^ i^e^^l decay [10] using a stopped pion beam. More than 
40,000 TT"*" — > z/ee"'"7 events have been collected, allowing for a very precise measurement 
of the branching ratio. In this experiment, a more significant discrepancy (about 8a [35]) 
between data and SM predictions has been reported in the kinematical region of high- 
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energy photon/low-energy positron. A significant number of expected events are missing. 
As for tlie ISTRA anomaly, agreement with data can be improved by adding a negative 
tensor term Ft — —0.002, a bit smaller (in magnitude) than the corresponding one in [14]. 
More detailed analysis about the PIBETA experiment can be found in [36] and references 
therein. 

Now we analyze the impact of a tensorial coupling on the photon polarization asym- 
metry dA^/dx in radiative pion and kaon decays. In order to simplify the analysis, we 
will assume an universal tensorial interaction in both AS* = 0, 1 processes, parametrizing 
all the effects in a phenomenological coupling Ft as follows 

Mt = ^^Ft K,ef [e(7^.(l - 75)i^e] , (58) 

where q = d, and q = s for pion and kaon decays respectively. As discussed in section 
1, the photon emitted by the tensorial amplitude Mt in Eq.(58) is purely left-handed. 
This property can also be checked by noticing that the interference between Mt and Msd 
terms is proportional to the V — A combination [14]. Below we provide the additional 
tensorial contributions to the (photon) polarized Dalitz plot density. In particular, the 
following term pj^^\x, A) should be added to p^~^\x, A) in Eq.(17) 

P^t'\x:X) = 2AsoFt (ftUt{x,\) + 2^il^ f,,,{x,\) + ^i{V - A)fsoT{x,X) \ (59) 

where [11] 

UrixA) = \x\l-\), /,,,(a;,A)=a;(l+n-A-^), 
fsurixA) = x^(l-A) . (60) 



In Fig. 16 we show the dA^/dx asymmetry versus for two representative values of 
Ft = ±10"^ and cut A > 0.3, for the 7r+ h'e^'^l (left plot) and t'eC+T (right 

plot). As we can seen from these results the shape of photon asymmetry is quite sensitive 
to a tensorial coupling in the range of \Ft\ — 10~^. In particular, in the pion case, 
this sensitivity is more pronounced, and the position of zeros of the asymmetry strongly 
depend on Ft- On the other hand, in the kaon decay, large deviations should appear 
only in the region of large x, where the tensorial effect is enhanced. This difference in 
the two decays can be explained due to the fact that the tensorial interference, which is 
the dominant effect when Ft < 10~^, is always chiral suppressed, being proportional to 
TTie/iTiM- Thus, in the radiative kaon decay, this effect is more suppressed than in the pion 
case, due to the larger meson mass. We have explicitly checked that, in the corresponding 
pion and kaon decays in muon channel, the sensitivity of the asymmetry to Ft is very 
modest and we do not show the corresponding results. In conclusion, we suggest that the 
possibility to measure the photon polarization in pion, or even in kaon decays, could be 
very useful to clarify the controversial question of tensorial 
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Figure 16: Photon polarization asymmetry dA^/dx versus x and with A > 0.3, for two values 
of tensorial coupling Ft = ±10"^. Left and right plots correspond to vr^ Vf,e^^ and — > 
V(,e^"f decays respectively. The dark dashed curves stand for the standard model case {Ft = 0). 

7 Cancellation of mass singularities 

In this section we discuss the mechanism of mass singularities cancellation and the way it 
takes place in meson and muon polarized radiative decays. As will be seen a new peculiar 
cancellation pattern shows up in the particular case of the polarized amplitudes differently 
from the well known cancellation taking place in the inclusive unpolarized amplitudes. 

In a theory with massless particles a crucial test of the consistency of the computation 
is represented by the absence of mass singularities in any obtained physical quantity. Mass 
singularities are of two types: infrared and coUinear. Infrared divergences originate from 
massless particles with a vanishing momentum in the small energy soft limit. Physical 
states as, for example, a single charged particle, are degenerate with states made by 
the same particle accompanied by soft photons. This corresponds to the impossibility 
of distinguishing a charged particle from the one accompanied by given number of soft 
photons due to the finite resolution of any experimental apparatus. An infrared divergence 
appears in QED when the energy of the photon goes to zero as a factor of the form: 

Jo e 

where e = -g^ is the fraction of the energy of the photon with respect to the total available 
energy E for the process. The Bloch Nordsiek theorem [37] assures the cancellation of 
infrared divergences in any inclusive cross section. CoUinear divergences, instead, come 
from massless particles having a vanishing value of the relative emission angle. In QED, 
specifically, when one or more photons, in the limit of zero fermion mass, are in a coUinear 
configuration i.e. with emission angle ^ 0. Physical states containing a massless 
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charged particle are degenerate with states containing the same particle and a number 
of coUinear photons. Any experimental apparatus, having a finite angular resolution, 
cannot distinguish between them. The angular separation of two massless particles with 
momenta p and k is such that they move parallel to each other with a combined invariant 
mass for ^ — > 0: 

q^^{p + kf = 2po E^{1 - cos ^) ^ (62) 

even though neither p nor are soft. Here 9 is the emission angle of a photon with 
respect to the fermion. The inclusive procedure of integrating over the photon emission 
angles by keeping the fermion mass finite does not give rise to any collinear singularity. 
The divergence appears in the limit ^ ^ as the presence of a logarithm of the form 
log(|) ^ log(^). 

For collinear singularities, as well as for infrared ones, the case for inclusive unpolar- 
ized processes is well known and it is governed by the Kinoshita-Lee-Nauenberg (KLN) 
theorem [38]. For the collinear singularities the KLN theorem guarantees that coUinear 
divergences cancel out if one performs a sum of the amphtude over all the sets of degen- 
erate states order by order in the perturbative expansion. For the amplitude of a single 
photon emission a combination of collinear and infrared singularities gives, for instance, 
contributions of the type: 

„ a de d9 , 

TT Jo e Jo 1 — COS B 

In order to discuss the above aspects on the cancellation of lepton mass singularity in 

the polarized pion decay, we first recall the mechanism taking place in the unpolarized 
case [39, 40]. In general, the decay rate is made free from mass singularities in the ordinary 
way: the cancellation of divergences occurs in the total inclusive decay rate at order 0{a), 
namely in the pion case 

r(incl) ^ ^ ^ ^ ^g^^ ^ ^g4^ 

when the full 0{a) order contributions are included, i.e. those relative to real and virtual 
photon emission [39, 40]. However, for a pointlike (structureless) pion, due to the chirality 
flip of flnal charged lepton, the pion decay amplitude is always proportional to m; and 
vanishes in the — > limit. In other words, in the limit — > the decay rate is made 
finite from mass singularities in a trivial way. For example, as we will see later on, a term 
proportional to Log(m/) will remain in the inclusive width due to the mass renormalization 
of the charged lepton in the virtual contributions to r(7r ve). However, since it will 
be multiplied by mf, it will give no troubles since Fq ^ tends to zero at the same 
time. However, as pointed out by Kinoshita in [39], the leading \og{mi) terms in the IB 
contribution to r(7r — > 1/67) cancel out exactly when one adds the virtual contributions. 
In other words, the mass singularities in the \og{rni) terms should cancel independently 
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from the fact that the effective couphng in the pion decay is proportional to the charged 
lepton mass or not. The cancellation mechanism of these Log terms shows a non trivial 
aspect of the KLN theorem in the pion decay. For this reason, in the following discussion 
we will consider the following ratios r(7r — > ve'yj/VQ and r(7r — > uej/VQ which survives 
the limit m; — > 0. 

Let us now consider the case of radiative polarized decays within the soft and collinear 
region for the radiated photon. We will investigate in this section the mechanism which 
will assure the finiteness of the lepton distribution against the appearance of infrared and 
collinear singularities on the above ratios of widths. Let us start by the inclusive distribu- 
tions in terms of the final lepton energy y as listed in Eq.(71). The Inner Bremsstrahlung 
contribution contained in Eq.(74) in the — > limit is composed by the four expressions 
corresponding to the various polarization states of the final photon and lepton respec- 
tively. The last, RR polarized term is identically zero. The remaining three are related 
to the left-handed ( first and third ) and right-handed ( second ) lepton respectively. 

The logarithms Li, L2 do correspond to collinear contributions. By integrating the 
double-inclusive distribution of Eq.(68) one gets in the expression of Eq.(72) for the IB 
case that the expressions for depend on the logarithms Li and L2 respectively 



The first logarithm represents the case of the photon being parallel to the lepton, the 
second collinear logarithm for mi ^ and ttim — ^ corresponds to the case where the 
photon is parallel to the decaying meson [3]. Clearly, it is only Li which is affected by 
the true collinear divergence in the limit — > 0. 

With respect to the unpolarized inclusive case some differences are worth to be noticed 
here: 

• Different polarization amplitudes do represent independent observables in the decay. 
Therefore if we consider the two cases of a right-handed and left-handed lepton they 
have to be also separately finite. 

• At zero order in the pion decay the angular momentum conservation imposes to the 
lepton to be left-handed. By radiating a photon a total zero angular momentum is 
assigned to the final state even if a right-handed lepton emits a left-handed polarized 
photon. This contribution is represented by the second term in Eq.(74). 




These terms do give rise to two kinds of "collinear" logarithms: 



L 



1 
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• For 2/^1 only the second term in Eq.(74), corresponding to the LR polarization, 
contribution is finite, i.e. it is zero: 

1 dvff^ a 



lim '-^ — = — (l — ii] = . 



n^o 2/^1 To dy 27r ^ 

This fact shows that the LR term, corresponding to the anomalous term [41], it is finite by 
itself without the need of any cancellation mechanism in the infrared y ^ 1 and collinear 
limit r; — > 0. A detailed discussion of the undergoing dynamics can be found in Ref.[3] 

Analogous conclusions, regarding the finiteness of the right-handed lepton contribution 
in the — ,|/ — > 1 limits to the structure dependent terms and the IB x SD, 

see Eqs.(75) and (76) respectively, hold there as well. 

Let us now consider the contributions of the type LL and RL giving rise to a left 
handed lepton. Manifestly the first and third term of Eq.(74) are divergent in the ri — > 
, y — > 1 limits. The expression obtained by adding first and third contributions in 
Eq.(74) is: 



1 



r„L dy +^J=^^[(l + !'-i.-i^) + (!'(!'+l)-i.!'' + i2(l-2,)) 

which is divergent both in the collinear and in the infrared limit. The coefficients of the 
collinear logarithms remain different from zero as r; — > and y — > 1, leaving to a divergent 
expression. As for the unpolarized case for the left-handed lepton contributions one needs 
to consider the additional virtual contributions in order to cancel infrared singularities 
[39]. 

The case of the left-handed lepton includes also the diagram of the virtual photon i.e. 
the one with a photon line connecting meson and charged lepton. This diagram docs not 
add any angular momentum to the zeroth order term since a virtual particle does not add 
angular momentum to the final state. The amplitude containing the virtual photon gives 
rise, therefore, to a lepton neutrino final state having the same helicities as the ones of 
the tree level amplitude. The combination of real and virtual contributions should, in the 
left-handed lepton channel, add among each other to give a finite result. This mechanism 
is the same taking place for the cancellation of singularities for the inclusive, unpolarized 
amplitudes as we will discuss in more details below. 

According to [39] the total width for the unpolarized IB contribution to r(7r — > ue^j) 
is given by 



^ib{xo) a 
— = -\b{ri} 



xq 1 3\ 

log ^ - iog(i - n) - - log(n) + -j 



n(10-7n), , , , 2(l + n),,, , , 15-21r, 



log(n) + ^L(l-..) + ^}, (65) 



A{l-rif ^""^"^ ' 1-ri ' 7(1 -n) 

where Xq is the minimum photon energy which regularizes the infrared divergence in the 

l+x 
(l-x) 



photon mass, the function b{x) = log(a;) -|- 2, and L{x) = /o^log(l — t)dt/t. For a 
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generalization of the result in Eq.(65) to the inclusion of the leading logarithmic terms to 
all orders in perturbation theory see Ref. [42-44]. 

For the virtual 1-loop contribution one has to consider the radiative corrections to the 
operator grriiQ, where Q = ■0/(1 — "y^jipufn, with (p^^ is the pion field and m° is the 'bare' 
mass of the charged lepton. These corrections split in two separate contributions: F^^) 
given by the correction to the operator Q and F^^-* arising when one try to express the 
bare mass m° in terms of the renormalized lepton mass mi, namely m° — mi — 5mi with 
5mi = Ifm/ (log(A/m/) + \) [40]. For F^^), one has [39] 

_ = _ (,og(^) - j,og(n) + j) + ^loglO + 5} + 5J (66) 

Notice that the last term, containing the ultraviolet cut-off A needed to regularize the UV 
divergency, can in principle be absorbed in a re-definition of /„■ at order a, see Ref. [40] 
for more details. 

As can be seen by comparing the results in Eqs.(65) and (66), the logr^ terms surviving 
the limit ri cancel out in the sum of virtual and real emission contributions as a 
consequence of the KNL theorem. Finally, for the total contribution to the unpolarized 
inclusive decay rate at order a, including the contribution of F^^-*, one gets [39,40]: 

F('"'=') ^ a r3, , , 13 TT^ . 

_ = l + -{-log(n) + ---}, (67) 

where we retained only the leading terms in me limit. As previously mentioned, the 
appearance of the logr; term in (67) is due to the renormalization of the charged lepton 
mass which does not follow the same pattern of coUinear mass singularities discussed 
above [39]. For simplicity, we omitted in (67) the term containing a log(A/m,r), since it 
can be absorbed into a re-definition of /„■ at order a inside Fq. 

As stated above, in the right-handed case, on the contrary, the mass singularities 
cancellation occurs with a different mechanism. Infrared and collinear limits in the ratio 
r/s^Vro give separately a finite result. In particular, the coefficient of the collinear 
logarithms for the right handed lepton case is the lepton mass, instead of the usual 
correction factor coming from the soft and the virtual photon contributions. 

The particular cancellation mechanism occurring in the right-handed radiative decay 
is originated by the combined constraints of the angular momentum conservation in the 
pion vertex and the one of the helicity flip in the photon-lepton vertex [3] . 

Let us now consider the case of the muon decay. As shown in Eq.(40) as for the 
meson case also in the muon decay lepton distribution we see that the Li?-photon-lepton 
polarized distribution is free from collinear and infrared singularities and goes to zero 
in the infrared limit y ^ 1. The remaining RL and LL distributions, apart from the 
identically zero RR term, do give a finite contributions in the y I limit, provided that 
the same xq cut-off is also set free to go to the soft kinematical limit i.e. xo — >■ 0. In the 
muon case the pattern of singularities cancellation repeats itself as for the meson case. 
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8 Conclusions 



We have computed polarized distributions in radiative meson and muon decays, by tak- 
ing into account final lepton and photon hclicity degrees of freedom. The definition of 
photon polarization asymmetry has been introduced, allowing a new approach to investi- 
gate interaction dynamics via a finite and universal quantity directly associated to parity 
violation. Analytical and numerical results for the polarized distributions and branching 
ratios, as well as for the photon polarization asymmetry, have been explicitly derived. 

The main results of the photon polarization analysis in meson decays, inclusive in the 
spin degrees of freedom of the final lepton, can be summarized as follows. In the pion 
case, the production of hard photons in association with soft positrons, are mainly favored 
to be left-handed polarized. However, when the positron energy increases, the relative 
gap between Icft-and right-handed photons decreases, due to the increasing contributions 
of hadronic structure dependent terms. Remarkably, in the kaon decay, when energy cuts 
E-y > 25 MeV and Ee+ > 120 MeV are imposed, both photon and positron are mainly 
right-handed polarized and a large and negative photon polarization asymmetry is ex- 
pected. Regarding the corresponding meson decays in muon channel, the left-handed 
photon production always gives the leading effect. The same behavior is observed in the 
radiative muon decay ii~ — > f^e^Pg 7. All these results can be easily explained in terms 
of angular momentum conservation and parity violation. 

We have also systematically analyzed the mechanisms of cancellation of infrared and 
coUinear divergences in polarized meson and muon decays. It has been shown that the 
finiteness of the polarized amplitudes takes place in a different way for left- with respect 
to right-handed final leptons when inclusive results in the photon polarization degrees of 
freedom are taken into account. 

Finally, we propose a possible test using photon polarization in order to solve the 
controversial issue of large tensorial couplings in lepton-quark interactions, as suggested 
by the recent observed anomaly at the PIBETA experiment. In particular, it is argued 
that the measurement of the photon polarization asymmetry may constitute a sensible 
test to resolve such controversial issue in radiative pion decay, providing a sensitive probe 
to hadronic form factors as well as to new physics effects in meson radiative decays. 

We believe that all these new results could open a more extended perspective into 
the physics of the scmilcptonic weak decays. In particular, the less inclusive approach to 
the polarized processes, by explicitly taking into account lepton and photon polarization 
degrees of freedom, could allow, when the experimental conditions make it compatible, 
a new quantitative approach and a more detailed inspection of meson and muon decays. 
Moreover, we are confident that Standard Model physics as well as signals of physics 
beyond the Standard Model could be put under scrutiny and more closely investigated 
by using tests involving polarized quantities. 
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Appendix A 



Here we provide the most general results for the polarized meson radiative decay in the 
meson rest frame, as a function of x^y variables defined in Eq.(ll), and for both charged 
lepton I — e, n and photon hclicitics A; and A^, respectively. For later convenience, we will 
use the same notation adopted in section 2, where the symbols L and R are associated 
to particle helicities A = — 1 and A = 1 respectively. In particular, for the differential 
radiative decay rate normalized to its non radiative one, we obtain 



a 



1 



1 (pr'^^i^^i') 



Fo dxdy 27r (1 - n)^ 



P 



x,y) 



(68) 



where 

P^'"^'\x,y) 
P^'''''\x,y) 



Jm 4 ri Jm 



/ir'H^, y) + "^^^^^n^Hx, y) + ^{v + A) ^')(x, y) 

JM 



where the functions fi'B'^'^\x,y), fsD'^'^\x,y), and fiNT''^'{x,y) are given by^ 



n'n''\x,y) 



f^V'\x,y) 



l-y + ri 



2Aix^z^ {A; (x + y-l + n i-3x-y + xy) + rf) 
— Xiy{x + y— l)+ri(^4 — 6x — 6y + 3xy + y^— X y^^ 
+ rf{i + 2x-y)]} 



1-y + ri 



{Ai {x - 1 + ri) - Xi[{x - 1) y - ri {2x + y - 4)]} 



2Ai x^ ^2 
X {{x-1) {x + y-l)+ri} 



{Ai {{x-1) {y-l) + ri) 



1-y + ri 

2Ai 

\[y {x + y-i-xy)+ri {2x-y)]} 

{x - l){x + y- ^)+ri ( . , . 
— [Ai{l-x-y + n) 

>^i[y {x + y - 1) - n {2x + y)]} 

i-=^^{A^ {l-2x-y + xy + ri) 



^The symbol A; = \/y^ — 4,r'f appearing below should not be confused with the axial form factor A. 
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- Xi 2-2x-3y + 2xy + y^ - xy^ + ri (2 + 2a;-|/)]} 



1-y + ri 
2 A, xz 



{Ai + Xi[2x + y-2]}{{l-x) {x + y-l)-ri} (69) 



where A; = 1 and A; = — 1 correspond to right- (R) and left-handed (L) fermion polariza- 
tion respectively, and symbols Fq = r(M+ — > uil'^), z — x+y — l — ri, and A; = y/y^ — Ari. 
By integrating equations above in the photon energy x range 



l-\{y + Ai) < X < l-i(y-AO 
< y <l+ri, 



we obtain 



1 dV^^'^'^''^ 



a 



dy 



27r (1 - n) 



(70) 



(71) 



where 

F(R'^i){y) 



(y) 



ml {V - Af 



Fg''^\y) + ^{V-A) 

JM ^ A 



2Ai{l-y + ri) fj, 4n 



{y) 



i'lB 



(y) 



+ 



2Ai{l-y + Ti) ' 4n ■ /m ■ 2Ai 



M 



INT 



{y) 



where the functions F/^'^'^'^(y) are given by 

F(^'^')(y) = A,{(L2 + L2) (l + r(l-y))+A; (I/-3 + 

+ A,{(Li + L2) (2r2-y + r (2 + y - y^)) + (2 + r (t/ - 6) - y + y^)} 

Ff'^'Hz/) = A,{Li (r + r2-3ry + y') -(L2 (l + r-2y + ry)) 

- A; (1 - 3r + y) } + A; {Li (2r2 + 2r^ + ry - 5r^y + 3ry^ - y^) 

+ L2 (y-2y2-r (2-7y + y2) -2r2) +A; (y + y^ - 6r + 2r' - ry) } 



^[{y-lfy-r (y + y2_2)+2r^}-h^{(y-l) ((y-1) y' 



- r (4y + y2-4) +4r^)} 



^^'Hy) - ^(y^-2ry(2 + y)+8r^) + ^(y^-4r) 



A; 



24 



24 



F^N'rHy) = (l-y + r){A,(L2 + Li(y-l))+A,(L2(y-2) + Li((y-l)y-2r))} 



47 



F^Nr'^'Ky) = -il-y + r){Ai{Ai + L2 + L,r) +Xi{2Ai{r-l) + L,r {2r-y) 

+ L,{y-2))} (72) 

where 

f y + Ai-2ri \ f y + Ai - 2 \ 

Now we expand the formulas above in the large electron energy region y ^ y/r. By 
retaining only the leading terms in ri expansion we obtain 

• I IBP contribution 



1 dr£'^^ a 



lim 

ri-*o To ay ztt y 

1 dr£'^) a / N 

hm ; = — [1 — y) 



^(l + y-Li-La) 



To 27r 

lim^^^ = ^^{y{y + l)-L,y' + U{l-2y)) 

lim -^^^ = (74) 
n-o To dy ^ ^ 



• |SD|^ contribution 

hm 



n-o r„ rfy 2^/1, 8r, " ^ 



1 dri^) _ a ml {V + Af , 

• IB X SD + c.c. contribution 

1 (iri!^T^ a ruM f l-y 

hm ; = — [V — A) L2 

n^o To dy 2% Jm y 

hm = -—-—{V-A)[L2 + Liy] 

n-^o To dy 2tt /m ^ ^ y 

J^of:^^ = -^(\/ + ^)(2L, + y(2 + y)) 



n-^o To 27r /m ^ ^ 2y 

n-»o To dy 27r /m ^ ^ 2y 
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where Li^2 = ^i^n^o (-^1,2) and so 



L2 = 



21ogy-log(l-y) -logn + 0(r) 
log(l-?/) + e»(r). 



(77) 



Notice that the Icpton mass inside Li is needed in order to regularize the colhnear diver- 
gences. The above results in Eqs.(74)-(76) are obtained in the approximation r/y^ <^ 1 
and are not valid near the region of minimum y ~ ^/r. Nevertheless, there is always a 
real infrared singularity in the photon energy spectrum, which is present in the terms 
log(l — y) when y — > 1 even if the electron mass is taken into account, corresponding to 
the known soft photon singularity. This divergent term for y — > 1 is necessary in order to 
cancel the infrared singularity appearing in the one-loop corrections of non-radiative de- 
cay, as required by the KLN theorem. Details of the cancellation mechanism for polarized 
decays are reported in section 7. 
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Appendix B 

In this appendix we provide the expressions for the basic functions go^i, ^0,1 ; (jio,1) ^-nd 
(jo,i appearing in section 2 for the differential decay rates of radiative muon decay. In the 
muon rest frame, the differential decay width is given in Eq.(37) 

A {90 + ^0) + Ae {gi + A^ Qi) (78) 



To dx dy dcosO Stt x 
where ;s = — Ag cos^), Ag = \Jy^ — 4r, and the functions 510,1 and ^0,1 are 

+ 2^(-3?/ + 2|/2 + 3^_4y^ + 2^2)+a;(^6y^-4y3-62 + 2|/z + 8y2 2 

- -6y^2 + 2^^)} +r{4a;^ + x2 (-6 + 87/-5;s) +2^2 (4_3^ + 3^) 

+ 2x^ (-3y + 2y^ + 3 z - y z + z"^) + X z (-Sy + 6y'^ - 6 z - 6y z + 3 z'^)} 

+ 2r^x^ (4:-3x-3y + 3z^ 

go = z[2x'^ + x^ {-3 + 6y-2z) + 2z^ {l-2y + 2z) 

+ (-6 y + Ay"^ - 2 z - 2 z'^^ + X z (g - 6y + Ay^ + z - 6y z + 2 z'^^} 

- r{4x'^ + x {-2 + 6y -3z) -6z^ -2x'^z {l-3y + z) 
+ x^ {-6 + Ay + z)} + 2r'^ x"^ {x + 3 z) 

gi = yz[2x'^ + x^ {-3 + 6y-2z) + 2x'^ (^-3y + Ay^-z-Ayz + z'^^ 

- 2z(^-3y + 2y^ + 3z-Ayz + 2z^) 

+ x (^-6y^ + Ay^ + 6 z - 2y z - Sy^ z + 5 z'^ + Qy z^ - 2 z^^ } 

+ r { - 4x^ + x^ {6 - 8y + 4 z) + 2 z"^ (^-12 + 8y + y^ - 8 z - y z) 

+ x^ (6y -8y^ - 4z + 7yz - 4z'^) +XZ (24y - 16y^ -2y^ + 16z + 18yz 

+ 2y'^z-yz^^+x'^ (dy'^ - 4y^ - 22y z + 6y'^ z + 6 z^ - lOy z'^ + 4z^^} 

+ 2r^{x'^ + 4xyz-4z'^ + x^ {8 + y + 2z) 

+ x"^ (^-12 + 8y + y^ -8z-yz + z'^)} - 
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y z {2x^ + x^ {-3 + 6y - 2 z) + 2 z^ {1 - 2y + 2 z) 
(^-6y + 4y^ - 2 z - 2 z"^) + X z (6 - 6y + 4y^ + z - 6y z + 2 z'^)} 
r[-4x^ + 2 {4 + y) z^ + x^ {6-8y + 4z) 
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+ xz^ (-8 - 2y - 2y'^ + y -2x^ z (^-5 y + 5y'^ + z-5yz + 2z'^^ 

+ (6y-Ay^ -Az-Syz + Az"^)} -2r^x'^[3x'^ + x {-A + 3y-2z) 

+ {A + y-z)z}. (79) 



After integrating the above distributions in cos^^, on the range —1 < cos^^ < 1, wc get 
1 dVib^^"^ a 1 



Go + Go + Ag {Gi + A-y Gi) 



To dxdy 247r AeX 

where the functions Go,i and Go,! are given by 

Go = -6A,L{2x + 2y-3) [x^ + 2xy + 2y^)+y^{2Ay {2y-3) 

+ 6x (A + y) {Ay-3)+x^ (36 + (33 - lOy) y) + 2x^ (6 + y (2y-3))} 

+ r{-6x (AeL (5x-6) + 8 (x (3 + x)-6))-12(-24 + AeL (4 + x) 

+ x (26 + (11 - 2 x) x) ) 1/ + 4 (9 Ae L - 24 + x ((11 -5x) x- 42)) y^ 

+ 9 (x - 4) (2 + a;) 2/^}+4r^|24 (3 y - 4) + a; (9 Ae L + a; (4a;-9y-4) 

+ 18(4 + 1/))} 

Go = x{6A^L {x + 2y) {2x + 2y -3)+y^ (l2 {6- X {3 + x)) 

- 3 (18 + X (2 X - 1)) y + 2 (6 + X (2 X - 5)) y^) + r (48 {x {3 + x)- 6) 
+ 12 (18 + x (2 X - 1)) y - 4 (x - 1) (5 X - 6) y^ + 9 (x - 2) y^ 

- 6AeL (x + 6y-2)) +4r^(9AeL + 4 {x - 3) (2 + x) - 9 {x - 2) y)} 

Gi = Ae{l2 (3 + r-2a;) (x^ + r (x^-S)) +4[-3a; (x-l) (6 + x) 

+ r (48 + X (18 + X (4 X - 13)) )] y - 3 ( - 24 + r (x - 4) (2 + x) 

+ X (26 + (11 - 2 x) x) ) y2 - 2 (24 + X (12 + X (2 X - 5))) y^} 

+ 6L{2r^ (-4 + x) (2x-y) + y (2x + 2y-3) (x^ + 2xy + 2y^) 

+ r(4 (x-1) x2 + (24 + x (7x- 22)) y + 2 (3x - 8) y^ - 2y3)} 

Gi = x{Ae[l2(8r + (r-3) (1 + r) x-2 (r-1) x=^) +4(3 (x-1) (6 + x) 

+ r (x-4) (3 + 4x))y-3 (-18 + r (x-2)+x-2x^) y^ 



(80) 



- 2 (6 + x (2x-5)) y' 



+ 6L 



2r^ (2 X - y - 4) - y (x + 2 y) (2 x + 2 y - 3) 



+ r (-4x^ + 10 (y- 1) y + x (4 + 3y))]}, 



(81) 
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where L = log (^^^x^ji s-iid r = ■m'l/m'j^ . The above results in Eq.(80) completely agree 
with the corresponding ones in Ref.[5] after summing over the photon helicities. 
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